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ABSTRACT 

The potential existence of a giant planet orbiting within a few AU of a stellar 
remnant has profound implications for both the survival and possible regeneration of 
planets during post-main sequence stellar evolution. This paper reports Hubble Space 
Telescope Fine Guidance Sensor and U.S. Naval Observatory relative astrometry of 
GD66, a white dwarf thought to harbor a giant planet between 2 and 3AU based 
on stellar pulsation arrival times. Combined with existing infrared data, the precision 
measurements here rule out all stellar-mass and brown dwarf companions, implying 
that only a planet remains plausible, if orbital motion is indeed the cause of the 
variations in pulsation timing. 
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• . 1 INTRODUCTION 



Observations of extrasolar planetary systems in the post- 
main sequence have implications for the survival of Earth 
and the terrestrial planets, and indirectly measure the ro- 
bustness of planet formation processes. Precision radial ve- 
locity surveys are discovering a growing number of giant 
exoplanets orbiting post -main sequence stars, where the 



hosts are either subgiant (jjohnson et aL 
20luj) or first ascent giant stars (Gettel et al.l 



Niedzielski et al.l I200S ; iDollinger et al 



201ll ; lBowler et all 
20121; 



20091 ; ISato et all 



20081 ; lLovis k Mavorl (20071 ). These systems have under 



gone relatively mild evolution compared to later stages, 
and th eir planet pr operties may only reflect formation pro- 
cesses (|Curriell2009n . However, it appears the effect of star- 
planet tides with increasing stellar radius is non-negligible at 
this e arly stage l|Llovdll201ll ; iHansenllioiol ; IVillaver fc Livid 
2009), and eventually becomes crucial IjNordhaus et al 



2010). 



Candidate planets orbiting post-red giant branch stars 
are observed with vary i ng degrees of observat i onal evidence 
ilCharpinet et al.l l201ll ; ISetiawan et all |2010| ; ISilvotti et ail 
2007), but these few detections currently lack the benefit 
of independent confirmation (e.g., transits) and statistics 
that have corroborated the conventional exoplanet popula- 
tion. Planets orbiting these more highly evolved stars and 



stellar remnants provide tests of primordial planet forma- 
tion, long-term star-planet and pla net-disk evolution, and 
second-generation plan et formation (|Wickramasinghe et al.l 
|2010| ; iMelis et alJlHoOflh . The latter scenario i s the only vi- 
able mechanism to produce the pulsar planets l|Hansen et al.l 
2009), and they remain the only post-main sequence plan- 
ets conf irmed by independent means (grav itational pertur- 
bations; IWolszczanlll994 iRasio et alj|l992h . 

Planetary systems around white dwarf s offer a glimpse 
into possible futures of the So l ar System llVeras fc Wvattl 
120121 ; iDuncan fc Lissauerl Il998l ; ISackmann et alJll993h and 
the opportunity to study the composition of planetary solids 
l|Zuckerman et alj|201ol ; iFarihi et a"ill2009l ; |jurall200gf ). Fur- 
thermore, and perhaps surprisingly, because white dwarfs 
outnumber A and F-type stars in the solar neighborhood, 
they may represent the majority of the nearest planetary 
systems formed at intermediate-mass stars(!). While there 
have been several ground- an d space-based searches for giant 
planets around white dwar f s dHogan et al. 2009; Farihi et al.l 
120081 ; iMullallv et ail 120071 ; iDebes et al.l 120051 '). to date the 
only published ca ndidate comes from variations in pulsation 
timing of GD66 l|Mullallv et al.ll200Sl) . 

This paper describes interferometric and astrometric 
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constraints on stellar and low-mass companions to GD 66 
(WD 0517+307). The motivation for the study is described 
in §2, where the observational data supporting a compan- 
ion are reviewed, and theoretical considerations for planet 
survival are explored. The Hubble Space Telescope (HST) 
Fine Guidance Sensor (FGS) and U.S. Naval Observatory 
(USNO) observations are described in §3, and companion 
mass limits are derived from the analysis of these data com- 
bined with prior studies. 



2 MOTIVATION 

There are two lines of reasoning that led to a search for 
stellar-mass companions around GD 66, and each are dis- 
cussed in turn below. The first is the continued, increasing 
trend in the observed pulsation arrival times, and the second 
is the issue of planet survival within a few AU during the 
post-main sequence evolutionary phases of an intermediate- 
mass stellar host. 



2.1 Observational Considerations 

iMullallv et alj (|2008l ) performed photometric monitoring of 
15 ZZ Ceti stars over the course of four years and identi- 
fied a low-amplitude, sinusoidal variation in the expected 
arrival times of pulsations in GD 66. A single turnover in 
the observed minus calculated (or O — C) diagram was ob- 
served around epoch 2005.3. Assuming the pulsation fre- 
quency is perfectly stable, and that deviations in transit 
time are caused by orbital motion, the magnitude of this 
variation is given by 



Table 1. USNO Observations of GD66 over 10.25 yr 



cM 



(1) 



where M and m are the mass of the star and companion 
respectively, a is the semimajor axis, and i is the orbital 
inclination. This can be rewritten as function of the orbital 
period p using Kepler's third law. 



I G(M + m)p 2 



1/3 



(2) 



Thus, for a given periodicity in the pulsation arrival times, 
the amplitude is a linear function of companion mass for 
m < M. 

A sinusoidal fit to the GD 66 ti ming data was foun d 
with parameters r = 3.8 s, p = 4.5 yr ijMullallv et al J 12008). 
thus implying a = 2.4 AU and ms in^ = 2.2Afj up for 
M = 0.64 M Q (|Bergeron et al.l 120041 ). This solution pre- 
dicted a turnover in t he O — C diagram in late 2007 that 
did not occur. Rather, IMullallv et~al] (2009) later reported 
that the arrival times continued to increase and a revised 
fit yielded r ~ 5 s, p — 5.7 yr, and hence a — 2.7 AU, 
msini = 2.4A'/j up . This newe r fit predicted a tu rnover in 
2008 that again did not occur l|Hermes et al.|[2~010h . 

Based on this trend of increasing r, it became possible 
that the observed signal could be due to a stellar compan- 
ion. In this case, the initial turnover discovered in the O — C 
diagram would be a relatively low probability event. Despite 
being somewhat unlikely, this alternative merits investiga- 
tion due to the profound implications of a planet orbiting 



Astrometry: 






16.86 ± 0.16 mas 


""rcf 


0.97 ±0.12 mas 


Tabs 


17.83 ± 0.20 mas 


[lot 


+55.3 ± 0.1 mas yr -1 


/<<> 


-120.3 ±0.1 mas yr" 1 


Vt an 


35.2±0.4kms~ 1 


Photometry: 




V 


15.56 ± 0.02 mag 


B-V 


+0.13 ± 0.02 mag 


V -I 


-0.05 ± 0.02 mag 


M v 


11.82 ± 0.02 mag 



Note. Photometry is on the Johnson-Cousins system and given in 
Vega magnitudes. Conservative errors for these mean magnitudes 
are 0.02 mag. 

a stellar remnant. Stellar mass solutions to Equation [2] be- 
come plausible for r > 10 s and i > 85° , but existing infrared 
data (see §3.3) constrain such companions to be degenerate: 
brown dwarf, white dwarf, or neutron star. Notably, a sec- 
ond white dwarf can remain hidden in optical spectroscopy 
if at least of comparable mass to the primary. An exce l- 
lent example of this is PG 0901+140 (|Farihi et al.ll2005l ). 
a 3'.' 6 DA5+DA6 binary that exhibit s an apparently-single 
DA5.5 spectrum (jLiebert et al.|[2005t ). For such a compan- 
ion mass to be viable, t he increasing O — C trend observed 
by (|Mullallv et al.|[2009h would have to continue for several 
years. Interestingly, orbital separations of a few to several 
AU have been found for five double white dwarf systems us- 
ing H ST FGS observations l|Subasavage et all l2009l : iNelanl 
l2007h . lending credibility to a binary scenario. 



2.2 Theoretical Considerations 

Because GD 66 is a carbon-oxygen core white dwarf, it has 
passed through both the first ascent (RGB) and as ymp- 
totic giant branch (AGB) phases. feergeron et alj (|2004l ) give 
spectroscopically derived stellar parameters for GD 66 of 
T eB = 11980 ± 200 K, log g(cms' 2 ) = 8.05 ± 0.05, imply- 
ing M = 0.64 ± 0.03 Mq. Using the average of three dif- 
ferent initial-to-fin a l mass relations dWilliams et alj 120091 : 
iKalirai et ail 120081 : iDobbie et all |2006h yields a possible 
range of main-sequence progenitor masses between 2.2 and 
2.6 Mq. As a basic line of reasoning that favors a stellar 
companion over a planet at a few AU, consider that the 
maximum AGB radius for this range of prog enitor stars is 

in the vicinity of 3 AU (|Villaver fc Livid 120071). 

However, as explored in detail bv lNordhaus et all (|2010l ) 
there are at least three processes that determine the orbital 
fate of low-mass companions to post-main sequence stars 
such as whit e dwarfs: 1) orbital expansion due to mass loss 
(jJeansj ll924j) , 2) tidal dissipation of orbital energy, and 3) 
destruction or survival in a common envelope phase. It is 
well known that very low-mass stellar and brown dwarf com- 
panions are capable of surviving withi n a common envelope 
l|Maxted et al.ll2006l ; [Farihi et al.ll2005l ) tha t results from un - 
stable mass transfer from a giant primary (|Paczvnskilll976l ). 
and these short- period systems are manifes t among white 
dwarf binaries l|Schreiber fc Gansickd 120031 : ISchultz et al.l 
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GD 66 Seasonal Averages 
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Figure 1. Residuals obtained after fitting proper motion and 
parallax to the USNO astrometric observations of GD66. Shown 
are the seasonal averages and the standard error of the mean. A 
total of 172 frames taken on 151 nights were used in the parallax 
solution, with 9 to 32 observations per season. The behavior of 
the residuals in Declination are also observed in a field star near 
to GD66 on the CCD, and are thus systematic. 



Figure 2. Companion mass sensitivity for primary mass 0.66 Mq 
and a 3cr threshold of 3 mas deviation due to circular orbital mo- 
tion for GD 66. Plotted as dotted lines are the astrometric ex- 
cursion as a function of companion mass for several benchmark 
orbital periods. The dashed line corresponds to a detectable com- 
panion mass of 50 Mj up at nearly all orbital periods between 4 
and 20 yr. For circular motion, the companion mass sensitivity 
is immune to orbital inclination, and thus the USNO astrometry 
rule out a wide range of low-mass, stellar and substellar compan- 
ions capable of producing the O — C variations. 



1996). The inward pull of tidal torques and common enve- 
lope evolution, together with the outward expansion of or- 
bits beyond the grasp of frictional and tidal forces, effectively 
create a depopulated region of intermediate orbital separa- 
tions (Nordha us et al.ll2010l ). This bimodal distribution of 
low-mass, unevolved companions to white dwarfs, barren 
from ~ 0.1 to ~ 10 AU, has been empirically c onfirmed us- 
ing high-re solution optical imaging with HST (|Farihi et al.l 
|2010| . 2006) and by s pectroscopic and photometric mo nitor- 
ing from the ground l|Nebot Gomez -Moran ct al. 201 if). 



Following the prescription of Nordhaus et al, (2010), 



the relevant semimajor axis boundaries for single planets at 
GD 66 were computed including mass loss, tidal forces, and 
common envelope evolution. The main-sequence progenitor 
masses span the range of values expected for this system 
(see §3.1), while companion masses include the msini = 
2.4 Mjup determined from the pulsation timing analysis, and 
a 7 Mj up upper limit es timated from infrare d photometry 
and modeling (see 53.3; iMullallv et al.|[2009l) . Included in 
these calculations are the largest, initial semimajor axis for 
planets directly engulfed by the AGB star, and smallest, fi- 
nal semimajor axis for planets that avoid being swallowed. 
Planets are destroyed in the first instance, while in the lat- 
ter case their orbits are strongly influenced by tides yet just 
avoid the giant envelope. It is found that for a wide range 
of possible tidal prescriptions, all companions that avoid en- 
gulfment end up in a ^ 3.6 AU orbits, and all but the most 
finely tuned initial conditions lead to final separations sev- 
eral AU larger. Therefore, an extant planet at 2 — 3AU 
around a white dwarf remnant of a M ^ 2.2 M© main- 
sequence star would require an unexpected evolutionary sce- 
nario (e.g., capture or re-formation). These predictions im- 
ply that a planet is rather unlikely, and largely prompted a 
search for stellar-mass objects capable of producing the sin u- 
soidal timing variations observed by (IMullallv et al,l F2008). 



3 OBSERVATIONS AND DATA 
3.1 USNO Astrometry 

Optical CCD observations of GD 66 spanning just over a 
decade were carried out as part of the USNO faint star par- 
allax program. Imaging photometry was collected for the 
purposes of correcting for differential color refraction and to 
measure an absolute trigonometric parallax. In practice, be- 
cause observations were taken within (and usually much less 
than) one hour of meridian, differential color refraction was 
minimal. A complete discussion of the astrometric data ac - 
quisition and reduction can be found in iDahn et"al (2002). 
Briefly, the observations were taken using the USNO 1.55 m 
Strand Astrometric Reflector and Tektronix 2048 x 2048 
camera with 24.0 /im pixels at Cf! 325 pixel -1 . These obser- 
vations e mployed a wide fl-b and filter, similar to that de- 
scribed in iMonet et alj |l992). Independent photometry in 
the Johnson BV and Cousins I band was collected on two 
nights at the USNO 40 in tele s cope. Photometric standards 
of various colors from lLandoItl (|l992l ) were taken at multiple 
airmasses to correct for extinction and color terms. 

The USNO astrometric results are listed in Table [TJ 
placing GD66 at d = 56.llg'gpc, and yielding an abso- 
lute visual magnitude of My = 11.82 mag. Comparing this 
with the spectroscopically derived value of My = 11.75 mag 
based on T cff = 11 980 K l|Bergeron et al.ll2004T ). the par- 
allax favors a slightly higher mass. Using this same ef- 
fective temperature, white dwarf atmospheric models yield 
log g = 8.09, M = .66 M , and a cooling age of 420 Myr 
l|Fontaine et al.l [200lh . Although the astrometric data pre- 
sented here do not constrain the temperature of GD 66, 
the the ZZ Ceti instabi lity strip is narrow at AT e g ~ 
1000 K for log g « 8.1 (|Gianninas et all 120051 ). To within 
the errors introduced by its few % photometric variabil- 
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Figure 3. Interference fringes for GD 66 along the FGSlr x- and y-axes. The fringes of the science target are shown as solid red lines and 
do not differ significantly from that of the calibration point source, BD+84 12, shown as dotted blue lines. At the bottom of each panel 
are shown the residuals in black, which arise from photometric noise, indicating that GD66 is unresolved. These data rule out stellar 
companions at Am = 2.44 mag in F538W with separations greater than 10 mas, or 0.56 AU at the trigonometric parallax distance. 



ity (|Fontaine et al.| [l985). the optical BVI photometry of 
GD 66 yields colors precisely as expected for a DA white 
dwarf of the published effective temperature. 

Figure Q] plots the astrometric residuals for each observ- 
ing season after subtracting the proper motion and parallax, 
demonstrating a few mas precision. This relative astrome- 
try is sensitive to orbital motion induced by an optically- 
dark companion. Notwithstanding the systematically noisier 
residuals in Declination relative to those in Right Ascension, 
the astrometric monitoring of GD 66 should have detected a 
3a deviation of 3 mas due to orbital motion in the plane of 
the sky. Figure [2] marks this detection threshold in a plot of 
the astrometric excursion as a function of companion mass, 
for several benchmark orbital periods in a circular binary 
configuration. These mass limits are independent of orbital 
inclination for circular orbits, and degrade only in the spe- 
cific case of an eccentric orbit at high inclination with line 
of apsides near to the line of sight (cj m 90°). 

Because the astrometric companion sensitivity only di- 
minishes for finely-tuned orbital parameters, the USNO data 
rule out all m > 50 Afj up (0.05 JWq), optically-faint compan- 
ions, including very low-mass stars, neutron stars and black 
holes, for orbital periods p > 4yr. Remarkably, a period- 
dependent range of relatively low, brown dwarf companion 
masses can be similarly ruled out for longer periods up to 
20 yr. However, these data provide little or no constraints on 
white dwarf companions in this mass range, as a binary of 
equal mass and brightness will exhibit no photocenter shift. 

3.2 FGS Interferometry 

GD 66 was observed during HST Cycle 18 on 6 October 2010 
by FGSlr in its high angular resolution Transfer mode, us- 
ing the F538W filter which covers 4500 - 7000A. In this 
mode FGSlr repeatedly scans an object and provides data 
from which interferen ce fringes alo ng its two orthogonal axes 
can be reconstructed i Nelan 21) 1 1 ). A relatively bright com- 
parison star known to be a point source at FGS resolu- 
tion, BD+84 12 was observed as a calibration source. Fig- 



ure [T] plots the FGS data for GD 66 and calibration star, 
and reveals the white dwarf is unresolved to 10 mas and 
Am = 2.44 mag. 

The FGS interferometry is sensitive to white dwarf com- 
panions at nearly all possible orbital separations and incli- 
nations. While two equal mass degenerates would be readily 
detected down to 0.56 AU and p — 0.4 yr, such short periods 
are not consistent with the p > 4yr, pulsation timing varia- 
tions of GD 66. For the periods fitted to the timing data of 
GD 66, Equation [1] dictates that a stellar-mass companion 
should be within several degrees of face-on. Therefore, the 
FGS data would have readily detected orbital separations 
from a few to several tens of AU, and thus rule out com- 
panions of comparable brightness and mass. However, in the 
very low probability event that the O — C turnover was a 1 in 
10 5 chance detection, the FGS data also constrain a putative 
double degenerate binary for any inclination. If the observed 
O — C minimum in 2005.3 corresponds to a binary system 
in conjunction as seen from Earth, then the projected sep- 
aration of the stars has been widening since, yet must still 
have an angular separation below 10 mas = 0.56 AU from its 
non-resolution by FGS in 2010.8. This implies a system pe- 
riod greater than 10 5 yr for edge-on orbits and total system 
masses M > 1 Mq . All double white dwarfs with shorter 
periods would have had a wider projected separation and 
been detected by the FGS. 

Importantly, white dwarf secondaries fainter than the 
primary by AU > 2.4 mag are not feasible by normal, un- 
perturbed stellar evolution; the total system age is insuffi- 
cient to achieve such low luminosities. From the 0.66 Mq 
white dwarf mass derived from the trigonometric paral- 
lax, the main-sequence progenitor of GD 66 had a mass be- 
tween 2.4 and 2.8 Mq (IWilliams et all 120091 ; IKalirai et all 
2008; Dobbie et al.l 120061 ), a hydrogen-burni ng lifetime in 
the range 460 - 660 Myr l|Hurlev et all l2000t ). and a total 
system age of 0.9 — 1.1 Gyr. Any cool and massive white 
dwarf (total age ~ cooling age) can become no fainter than 
M v = 13.9 mag over 1.1 Gyr (|Fontaine et al.1 120011 ) and 
would be detectable as AU < 2.1 mag. 
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3.3 Infrared Photometry 

Existing Spitzer IRAC photometry of GD 66 are sensi- 
tive to substellar companions that would n ot be detected 
by th e USNO or FGS optical astrometry. iMullallv et all 
(2009) used these infrared observations to estimate an up- 
per limit on spatially-unresolved (within 2'.'4 « 130 AU), 
self-luminous objects of 7Mj up , thus indicating the putative 
companion should have a planet-sized mass. Their method 
utilized observations of three white dwarfs with similar tem- 
peratures to constrain the photospheric flux ratio between 
3.6 and 4.5 /im, reporting a la total uncertainty of 0.6% in 
the ratio observed for GD 66. As an e xample of an alter- 
native method using IRAC photometry, iFarihi et al.l (|2008l ) 
establish substellar companion mass limits at 15% above the 
predicted (or measured) 4.5 /im photospheric flux. Applying 
this procedure to GD 66 with the USNO parallax distance, a 
total system age of 1 Gyr (see §3.2) , and assuming its 4.5 /im 
flux measurement |Mullallv et al.l l2009) is due to the photo- 
spher e, one obtains an upper limit of 11 Mj up l)Baraffe et al.l 
2003). If one insists on only a 10% excess at this wavelength, 
the upper limit drops to and 9 Mj up . It should be mentioned 
that all such analyses depend on atmospheric models, none 
of which have empirical constraints at these masses, as well 
as the adopted total system age. Regardless of the method, 
the infrared data rule out all but planetary masses, accord- 
ing to models. 



4 CONCLUSIONS 

The combined observational data on GD 66 limit any binary 
companions orbiting within several AU to planetary masses. 
Specifically, any long-term, increasing trend in the pulsa- 
tion arrival times cannot be due to stellar-mass secondaries, 
which include low-mass stars, white dwarfs, neutron stars, 
and black holes with periods longer than 4yr. The USNO 
relative astrometric monitoring of just over a decade rules 
out stellar-mass, dark companions with periods between 4 
and 20 yr, while the FGS observations rule out virtually all 
white dwarf companions, regardless of orbital inclination. 
Based on substellar cooling models, infrared data further 
restrict low-mass companions within a hundred AU to have 
planetary masses. It is noteworthy that the trigonometric 
parallax and infrared photometry by themselves rule out 
a range of double degenerate scenarios, but the astromet- 
ric monitoring and interferometry provide significantly more 
stringent limits to such binaries. 

If the observed timing data at GD 66 are due to orbital 
motion, these new and exiting data rule out a vast range of 
realistic companion masses, and strengthen the likelihood 
of a planet-sized mass as the cause. Because stellar evolu- 
tion models and resulting star-planet interactions indicate 
a planet within a few AU of an intermediate-mass star is 
not likely to survive to the white dwarf stage, the putative 
planets at GD 66 and other post-RGB stars V391 Peg and 
KPD 1943+1405 (KOI 55), if real, may have been dynami- 
cally injected or formed in a second-generation of planet for- 
mation. Further study of these and future, post-RGB planet 
candidates is needed to better understand the population 
and architecture of planetary systems around white dwarfs. 
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